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ABSTRACT 

We find that all classes of galaxies, ranging from disks to spheroids and from dwarf spheroidals to 
brightest cluster galaxies, lie on a two dimensional surface within the space defined by the logarithms of 
the half-light radius, r^, mean surface brightness within r^, le, and internal velocity, V^ = {\v1 +cr^), 
where Vc is the rotational velocity and a is the velocity dispersion. If these quantities are expressed 
in terms of kpc, L© pc~^, and km s~^, then this surface is described by the equation log re — logl^^ + 
log/e + logTg + O.S = 0, where we provide a fitting function for Te, the mass-to- light ratio within r^ in 
units of Mq/Lq, that depends only on V and If.. The scatter about this surface for our heterogeneous 
sample of 1925 galaxies is small (< 0.1 dex), and both the scatter within one of the galaxy subsamples 
(1319 disks) and the analysis of subsamples with independently derived mass-to-light ratios suggest 
that the intrinsic scatter could be as low as ~ 0.05 dex, or 10%, prior to any correction for observational 
errors. This small scatter has three possible implications for how gross galactic structure is affected 
by internal factors, such as stellar orbital structure, nuclear activity, or mass loss history, and by 
external factors, such as environment or accretion history. These factors either 1) play no role beyond 
generating some of the observed scatter, 2) move galaxies along the surface, or 3) balance each other 
to maintain this surface as the locus of galactic structure equilibria. We cast the behavior of Te in 
terms of the fraction of baryons converted to stars, 77, and the concentration of those stars within 
the dark matter halo, ^ = i?20o/'"e, where i?200 is the standard estimate of the virial radius. We 
derive expressions for 77 and ^, use an independent measurement of 77 to evaluate leading constant 
terms, obtain ry = 1.9 x 10~^(L/L*)T*F~^ and ^ = lAVr~-^, and relate these to each other via 
log?7 -|- log^ = — log Te + log T, + const. Finally, we present the distributions of 77 and ^ for the full 
range of galaxies and conclude that the high Te's of dSphs are due to low 77 rather than any differences 
in ^, that 77 is similar for spheroids and disks of a given 1/, and that 77 decreases with increasing V for 
systems with V > 30 km sec~^. For systems with internal velocities comparable to that of the Milky 
Way (149 <V< 163 km s^^) , rj — 0.14 ± 0.05, and f is, on average, ~ 5 times greater for spheroids 
than for disks. 
Subject headings: galaxies: formation — galaxies: structure 



1. INTRODUCTION 

Are galaxies fundamentally a simple family of collapsed 
objects, whose gross structure is describable by a few 
basic parameters, or are they highly complex systems 
whose structural properties are determined by a myriad 
of internal and external factors? 

If the former, there must be an analogous construct to 
the stellar Hertzprung-Russell diagram that testifies to 
deep, systemic structural patterns among galaxies and 
serves as a guide to a simple, if not entirely complete, 
analytic description of galactic structure. The study of 
stellar structure offers a beautiful example of the power 
of reductionism in astrophysics. By focusing on the HR 
diagram, investigators solved the problem of stellar struc- 
ture without needing to address other unsolved problems, 
such as the origin of the initial mass function. The ob- 
servation that the position of main sequence stars on the 
HR diagram is insensitive to their location in the Galaxy 
indicates that their structure does not depend sensitively 
on parameters that vary from one place to another. We 
now know that mass is the primary determinant of where 
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a star lies on the main sequence. Other physical char- 
acteristics, such as age, metallicity, and rotation, affect 
stellar colors and magnitudes (and therefore should be 
included in a complete model of stellar structure), but 
they are relatively minor factors along the main sequence 
and can be neglected in the interest of isolating the basic 
physics. 

Among galaxies there are hints of analogous "se- 
quences". These are referred to as ga laxy scaling laws 
and include the Faber- Jackson (FJ; [Fabcr & Jacksq 
19761) and TuUy-Fisher relations ( TF; iTuUv fc Fishe 
the Fundamental P lane (FP; iDjorgovski fc David 
Dressier et al .11 19871). an d the Fundamental Mani- 
a.b). Although it is not yet 
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evident that any of these is as fundamental for galaxies 
as the main sequence is for stars, they do imply that 
a limited number of parameters characterize the gross 
properties of at least certain subsets of galaxies. 

There are two arguments against using existing galaxy 
scaling laws as guides to a fuller description of galactic 
structure. First, existing scaling laws work only over a 
limited range of galaxy types and luminosities. While 
this is not an insurmountable obstacle — not all stars 
lie on the main sequence — it suggests that the cur- 
rent scaling laws are incomplete and that they will not 
astro'fi&4d\P ^ description of all galaxies. Second, for histor- 
leal reasons related to the empirical nature of the seal- 
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ing laws, their current formulation is not optimal with 
respect to possible theoretical constructs. For exam- 
ple, one determination of the F undamental Plane ha s 
logre = 1.24 log cr-0.82 log /e+7 (|J0rgensen et al.lll996D . 
where rg is the half-light radius, a is the velocity disper- 
sion, /e is the surface brightness within Tg, and 7 is a 
constant. It is unlikely that a simple theory would re- 
produce the 1.24 and 0.82 coefficients. 

In this paper, we attempt to address both of these 
shortcomings. Although we are not the ffist to hope to 
identify a unifying d escription of galaxies (see K-space; 
iBurstein et al.lTl997h . we achieve three of our key goals: 
1) to find an empirical relationship for all galaxies that 
has comparable scatter to those relations identified pre- 
viously for limited subsets of galaxies (TF, FJ, FP, and 
FM), 2) to isolate the critical additional knowledge be- 
yond the virial theorem that is needed to derive this rela- 
tionship, and 3) to begin constructing the bridge between 
a purely empirical relationship that utilizes observables 
and a theoretical one that is based on physical parame- 
ters. 

In summary, we begin from basic dynamical principles 
and examine the dimensionality of the family of galax- 
ies, ranging from dSph's to brightest cluster galaxies and 
from disks to spheroids. We address our basic question — 
how uniform are the gross structural properties of galax- 
ies? — by determining that a single scaling relation ex- 
ists that spans all luminosities and galaxy types and by 
quantifying its scatter. We employ an extended version 
of the fundamental plane formalism that reproduces the 
structural properties of all galaxies at a level compara- 
ble to that achieved with either the TF relation for disks 
or the FP for spheroids. We establish that connecting 
this new scaling relation to the virial theorem requires 
knowledge only of the mass-to-light ratio within r^, Tg, 
and that Tg can be accurately modeled as a function of 
the observed structural parameters themselves. We pro- 
ceed to describe Tg as a combination of the fraction of 
baryons converted into stars, 77, and the degree to which 
those stars are packed within the dark matter halo, ^. 
Using our empirical findings, we then calculate these two 
physical parameters for all of our galaxies and compare 77 
with independent measurements. The aim of this work is 
to define simple expressions for basic physical parameters 
of galaxies that may illustrate which physical processes 
drive the observed patterns of galactic structure, with 
the expectation that this will focus subsequent, more de- 
tailed theoretical work. 

2. THE DATA 

To determine the degree to which all galaxies are struc- 
turally similar, we need structural parameter measure- 
ments for galaxies ranging from spheroids to disks and 
giants to dwarfs. Part of the legacy of distinct scaling 
relations for different classes of galaxies are studies that 
provide the relevant information only for those partic- 
ular classes of galaxies. For example, there are exten- 
sive s tudies of spheroidal galaxies (e.g.. Ij0rgensen et al.l 
19961) that are entirel y distinct from those of spirals (e.g., 
Springob et al.l [20071 ). This dichotomy is partly due to 
the techniques necessary to measure the internal dynam- 
ics for disks and spheroids, but it also leads to the use 
of different photometric systems and definitions. It is 
impossible to resolve all of those differences, and many 



existing galaxy samples cannot be included here because 
they lack some necessary measurements. We describe the 
spheroid and disk samples that we use below. These con- 
stitute a heterogeneous dataset, but span the full range 
of galaxy types and luminosities, and require minimal 
corrections for internal comparisons. It is a testament to 
the robustness of our results that the many differences 
among the samples that we either ignore or only crudely 
correct (such as correcting the photometry to / band on 
the basis of average colors for different galaxy popula- 
tions) do not derail this investigation. 

2.1. Spheroids 

Since [Zaritskv et all (|2006bf ). there has been one key 
improvement in the av ailable data on low-mass spheroids. 
ISimon fc Gehal (|2007f ) present velocity dispersions, and 
a uniform set of structural parameters, for eight addi- 
tional Local Group dSph's, including some of the lowest 
luminosity sys t ems k nown. Adding these data to the 
IZaritskv et al.l ()2006a[ ) compilation greatly increases our 
sample for extreme values of luminosity, internal velocity, 
and effective radius. The lack of such data earlier pre- 
cluded our use of this range in the Gtting of the FM, and 
instead we showed that an extrapolation o f the FM accu- 
rately fit galaxies in this parameter range (| Zaritskv et al.l 
l2006b( l. Here, we fit to both th e previous da ta for the 
entire range of spheroid masses (jZaritskv et a l. 2006a, I3, 
and refere nces therein) and the new data for low mass 
spheroids (jSimon fc Gehal 120071 ). 

2.2. Disks 

We focus on t hre e particular disk samples: 

Pizaenoetall (l2007l) . ISpringob et"an (|2007l ). and 



Geha et all (|2006t ). Here we briefly describe the various 



data sets. 

Of the three samples, the lPizagno et all (|2007f ) sample 
allows the simplest comparison to the spheroid samples. 
The authors provide half-light radii, i— band magnitudes, 
and a range of velocity measures from their optical rota- 
tion curves. As they did for their TuUy-Fisher analysis, 
we use their Vso measurement, which is a measure of 
the rotation velocity at a radius that encloses 80% of 
the galaxy light. We correct the z— band magnitudes to 
Johnson by subtracting 0.4 mag (iFukugi ta et al.l[r995| ). 

The next simplest sample for comparison is that of 
fSpringob et al.l (|2007f ). who provide HI measurements 
of the rotation and /-band photometry. They do not 
tabulate half-light radii, so we calculate them based 
on the measures they do provide, the radius that en- 
closes 83% of the light and the radius of the 23.5 mag 
(sq. arcsec)"^ isophote, assuming an exponential surface 
brightness profile. Among galaxies for which all of the 
relevant data exist, we only reject systems with cz < 2500 
km sec~^, to avoid the loc al flow field. 

Lastly, the I Geha et al.l ([2006) sample is distinct be- 
cause it is primarily composed of low luminosity systems 
with very large gas mass frac tions. Because ga.s frac - 
tions are low in normal spirals (jRead fc Trenthamll2005| ). 
the gas can be ignored with little impact when study- 
ing scaling relations like TF for such spirals. However, 
studies of low mass galaxies show that accounting for all 
the baryons is critical in maintaining the s caling rela- 
tion (|McGaughet al.l 120001 : lMcGaughll2005t iGeha et all 
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I2006D . Therefore, we discuss the lGeha et"all ()2006D sam- 
ple separately in §3.41 We use their inclination- and 
turbulence-corrected velocities, and transform from r to 
/ magnitudes using the colors of late type spirals and the 
tabulations of iFukugita et alJ (|1995| ). 

3. RESULTS AND DISCUSSION 
3.1. Proceeding from the Virial Theorem 

In this section, we revisit the standard derivation of 
the FP to provide a framework and physical intuition 
for our observational results. We begin with the tensor 
virial theorem, use simplifying assumptions to rewrite 
the virial theorem in terms of observed quantities where 
possible, discuss the resulting equation and its implica- 
tions for the nature of galactic structure, and finally sug- 
gest a way to proceed even though some terms in the 
resulting equation cannot be expressed in terms of ob- 
served quantities. We apply this suggestion and explore 
it in quantitative detail in §3.21 and ^ 33. 31 

We begin with the tensor virial theorem, which is 
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where T and II are the contributions to the kinetic en- 
ergy tensor from the ordered and random motions, W is 
the potential energy tensor, and I is the moment of iner- 
tia tensor. In steady state, the left hand side of Eq. [l]is 
zero. We evaluate the trace of this equation and express 
the ordered component of the kinetic energy as ^AqMvI 

and the random as AiMa'^, where Vc is the circular ve- 
locity in disk galaxies, a is the line-of-sight velocity dis- 
persion for spheroidals, M is the mass of the system, and 
the A's represent the correction factors obtained by fully 
evaluating the appropriate integrals. Similarly, the po- 
tential energy is expressed as —BqGM'^/R, where i? is a 
characteristic radius that we define to be the half-light 
radius, r,,, and Bq is a correction factor obtained by fully 
evaluating the appropriate integral. Hence, without loss 
of generality, 



AqvI + Aia^ = Bo 



GM 



(2) 



All of the possible real- world complications are encapsu- 
lated in the yet unspecified j4's and in Bq, and, in prin- 
ciple, these could be extremely complicated functions of 
the formation history and environment of galaxies. The 
only assumption that we have made so far is that the 
virial theorem holds over these radii, which is reasonable 
for galaxies because re <C r^ir. 

To numerically evaluate Eq. [21 we now introduce two 
sets of simplifying assumptions that we will eventually 
test by determining whether we reproduce the observa- 
tions. First, is the kinematic simplification. We reduce 
the number of A parameters, by requiring Aov'^ + Aia'^ = 
AV^, where V = {^v'^ + a"^). This simplification is 
accurate if we are dealing with isothermal spheres and 
isotropic velocity dispersions. In such systems, at large 
radii, Vc for a purely rotationally-supported population 
equals ^/2a for a purely pressure-supported population. 
Furthermore, in such a pressure-supported system, the 
internal velocity dispersion is equal to the line-of-sight 
velocity dispersion. To evaluate V for disk galaxies, 
we will use the measured Vc and set a ~ 0, while for 
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Fig. 1. — The distribution of galaxies (444 spheroids, 1481 disks) 
in the (re,Ie,V) space. This projection was chosen by eye to 
demonstrate that galaxies lie on a fairly well defined surface. The 
data are discussed in [[2] The axes in the lower right show the 
orientation of the 3-space. The poorly populated tail consists of 
Local Group dwarfs and would be, we expect, well populated in a 
volume limited sample. Units for the three axes are kpc, Lq /pc^, 
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for disks and V = a for spheroids. 



spheroidal galaxies we will use the measured a and set 

Although the kinematic simplification relies on highly 
specific assumptions, both disks and spheroids satisfy the 
relevant conditions well, and this conversion has been 
used previou sly in various contexts (jBurstein et al.|[l997t 
iKassin et al.l[2007[ ). Optical disks are characterized by 
flat rotation curves, which imply that the mass profile is 
that of an isothermal sphere over these radii, and that the 
velocity tracers, H II regions or neutral hydr ogen, are on 
circular orbits (see lFaber fc Gallagheil[l979h . Spheroids 
also lie in mass distribut ions that are consist ent with be- 
ing isothermal spheres (jGavazzi et al.l l2007f ) . and their 
stellar velocity dispersi ons are nearly isotropic if the sys- 
tem is a slow rotator (jCappellari et al.ll2007f l. Due to 
the nature of our spheroid samples, we expect that there 
are strong selection biases against fast rotators among 
the more luminous systems (i.e., they would often be 
removed from Fundamental Plane studies) , and the low- 
est lu minosity systems show little rotation (jWalker et al.l 
l2007t l. Therefore, our sample is likely to satisfy the as- 
sumptions involved in the kinematic simplification, but 
we discuss possible signatures of failure in i j3.5l 

Continuing in our attempt to convert Eq. [2] into an 
equation that we can numerically evaluate, our second 
simplification involves the replacement of the difficult- 
to- measure M with Mg, the mass enclosed at r^- We 
refer to this as the mass simplification. We rewrite M^ as 
Tgie so that it is expressed as a function of observable 
quantities: L^ and Tg, the luminosity (= nr'^If,) and 
the mass-to-light ratio within Tg, respectively. We then 
replace Bq with B to account for the unknown difference 
between M and Afg. 

We now apply the kinematic and mass simplifications 
to Eq. [21 rewriting it as 



AV" = BG-KTerJe 



(3) 
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Finally, we take the logarithm of both sides and rearrange 
terms to obtain 

log Tg — log V^ + log /e + log Te — log A + log B = const. 

This equation leads to the rather dispiriting conclusion 
that galaxies populate at least a six dimensional param- 
eter space — more if yet unspecified parameters, such 
as Tg, are actually functions of additional parameters, 
like age, metallicity, formation history, bulge-to-disk ra- 
tio, or environment. Surprisingly, as shown in Figure [H 
galaxies populate a limited region of the (rg, /g, V) space, 
indicating a much lower dimensionality. 

One way in which the dimensionality of the galaxy 
family might be reduced from that suggested by Eq. S] 
is if Te, A, and B are functions only of re,Ie, and V. 
A simple variant of this scaling is referred to as "ho- 
mology" , in which the functional forms are assumed to 
be power laws. Because of the logarithms, the end ef- 
fect of rewriting Eq. [4] in such a variant is a change in 
the coefficients of the logT^^ and logJe terms. There- 
fore, the assumption of homology results in a prediction 
that galaxies lie on a plane in the (re,/e, V) space. The 
values of the coefficients describe the tilt of that plane. 
The success of the Fundamental Plane de s cription for gi- 
ant ellipticals (iDiorgovski fc DaviJIlQSTt iDressler et al.l 



HQSTtlJorgensen et al.lll996l : lBernardi et al.il2003l ) demon 
strates that, over the limited mass range of these galax- 
ies, the homology assumption holds surprisingly well. 
T his success was exten ded in the k— space formalism 
of iBurstein et al.l (|1997| ). where different classes of ob- 
jects were found to lie on different planes. However, 
the failure of a single plane to describe the distribu- 
tion of all spheroidal galaxies demonstrates that over 
a more extended mass range, which includes the most 
and least massive sp heroids, homology does not hold 
([Zaritskv et al.ll2006aD . 

Here we take a different approach in that we (1) assume 
that a77 galaxies (faint, luminous, disk, and spheroid) fall 
on a single manifold in the (re,/e, V) space and (2) ex- 
amine the behavior of Te that would make that possible. 
This approach is motivated by our earlier finding that the 
behavior of Te is bot h simple and qu alitatively reason- 
able for all spheroids (jZaritskv et al.l l2006a,b) — leading 
to what we termed the Fundamental Manifold (FM) of 
spheroids. So emboldened, we now assert that for all 
galaxies deviations from homology are dominated by the 
behavior of Te, ignoring variations in A and B among 
galaxies^. Our approach here represents a philosophi- 
cal departure from ours and others' earlier work, which 
usually focused on establishing or quantifying the tight 
empirically-derived scaling relationships (e.g., the FP or 
FM), because we posit the existence of a fundamental 
manifold of all galaxies and then examine the implica- 
tions. 

3.2. The Simplicity of Galaxies 

The treatment described in §3.11 and culminating in 
Eq. 0] is incomplete. The simple theoretical approach 
fails because it does not predict the low dimensionality of 

^ An alternate treatment that eventually leads to the same con- 
clusion is to group together Te, A, and B into a ge neric unknown, 
A, fit for A, and then use the argument in il3.3l to demonstrate 
that A oc Te. 



the data seen in Figure [T] On the other hand, the purely 
empirical treatment of fitting a manifold to the data in 
the (re, /e, V) space fails because it does not connect the 
actual functional form to a physical framework. Much 
like the case with the FP coefficients, directly fitting the 
data will subsume the behavior of Te, A, and B in Eq. [4] 
into the coefficients of the various structural terms (see 
axes in Figure [T]). Instead, we merge the two approaches 
by retaining the values of the coefficients derived from 
the virial theorem treatment as given in Eq. 31 assert 
that the most distinct break from homology occurs in 
Te, set A and B to be constants, and then solve for Te, 



log Te = log V — log /e — log Tg -l- const. 



(5) 



This approach may seem like only mathematical sleight- 
of-hand, but we will quantitatively test our association 
of Te with the dominant departures from homology in 



To proceed, we evaluate log Te — const, using Eq. [5] 
and plot the results in Figure [21 We then fit for the 
function, logT{, that describes these data and also plot 
the calculated values using this fit in Figure [2l Because 
we are fitting to a distribution of points in a 3-space, the 
fitting function will depend on two variables, and the nat- 
ural choices are those that are distance independent, V 
and If,- We want to minimize the fitting order, while still 
capturing t he behavior of t he dis tribution. As demon- 
strated by IZaritskv et al.l ()2006aD . there is at least a 
second-order dependence on log a and some dependence 
on log le , and so we fit to second-order in both log V and 
log/e and include cross-terms. Fo r this fit, we use only a 
randomly selected one-sixth of the lSpringob et all ()2007D 
sample to avoid having that sample dominate the fit. We 
present the coefficients of our fit in Table 1, but, because 
of the heterogeneous nature of the data and ou r avoid- 
ance of any type of Malmquist-like corrections (|Willickl 
Il994[ ). these numbers are far from definitive. 

The distinction between this work, with its complex 
characterization of Te, and either FP or K-space, with 
their assumption of homology, becomes evident when ex- 
amining Figure |2l The projections of the data and the 
fitting function in Figure |21 illustrate how even in projec- 
tion the functional form that describes Te deviates from 
power laws. The upper panels contain the inferred values 
of logTe — const, from Eq. jSlvs. either logF or log/e. 
The lower panels, which show the fitted values, logT{, 
and therefore have no intrinsic scatter, illustrate how the 
bulk of the observed scatter in the upper panels comes 
simply from the projection of a complicated surface onto 
these axes. In other words, the reason why galaxies of 
the same V have a range of Te's is not primarily because 
there is intrinsic scatter — say, due to age or metallicity 
— but rather because galaxies have a range of le- For a 
given V and /e, the scatter in Te is much smaller than 
that observed in Figure [21 To be specific, the scatter for 
the entire sample about the fit is 0.094 dex (24% rms in 
Te). In contrast, the observed scatter in logTe in the 
upper left panel of Figure |2l for 1.9 < y < 2 is 0.22 dex 
(66% rms in Te). 

We are now ready to evaluate the degree to which Eq. 
m describes our set of galaxies. Replacing logTe — const 
with logT{, we evaluate Eq. 5 and plot a rearrange- 
ment of the terms in Figure [31 By construction, Eq. |5l 
is satisfied on average when log Te — const, is replaced 
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Fig. 2. — Projections of logTe and logT^. We plot the projections of logTe — const., determined from Eq. \5\ vs. log 1/ and log/e 
in the upper panels. In the lower panels, we plot the values of logX^ for every galaxy in our sample using the fit given in Table 1. The 
lower panels illustrate how even with no intrinsic scatter in logX^ the projections show significant apparent scatter. We conclude that 
the bulk of the apparent scatter in the upper panels is due to the effects of projecting the complicated surface onto these axes rather than 
observational errors or intrinsic scatter. 

dence of Tg on accretion history, age, varying degrees of 
mass loss, or many other possible physical effects. These 
differences in Tg are not accounted for in T{, thereby 
potentially leading to a large scatter about the mean. 
However, the scatter is only 0.094 for the entire sam- 
ple. For reference, the scatter in this new relation for 
all galaxies is comparable to the scatter observed in ei- 
ther FP or TF studies for the relevant subset of galaxies. 
The scatter can be reduced slightly, to 0.084, if we cor- 
rect each galaxy sample separately for zero point shifts. 
These inferred zero point shifts, obtained by calculating 
the mean offset relative to the 1:1 line, are all comparable 
to plausible photometric errors, and generally correspond 
to a few hundredths of a magnitude. 

The success of placing all galaxies onto a single sur- 
face in the (rg, /e, V) space demonstrates that, to within 
the scatter (< 25%), the family of galaxies is a two pa- 
rameter sequence, i.e., measuring two of these structural 
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by log T{ , which is evident in Figure [3] The actual test 
of our approach comes from examining the scatter about 
the mean and whether distinct galaxy populations fall off 
the mean trend. If galactic structure depends strongly on 
parameters not included in this simple description, then 
the scatter will be large. In other words, two galaxies 
that are identical in the quantities V, r,,, and I,, could, in 
principle, have very different values of Te due to a depen- 
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Fig. 3. — Testing an equation of galactic structure, Eq. \E\ We replace logTe — const, with Ti to evaluate Eq. |5]for the entire sample 
(left panel) and for spheroids and disks separately (right panels). By construction, the mean relation should lie along the 1:1 line. The low 
scatter and the lack of systematic deviations for galaxy subsamples testifies to the universal nature of this simple relation. 



parameters specifies the third. This imphes that poten- 
tially important factors in galaxy development, such as 
environment, nuclear activity, star formation history, and 
accretion history, do not play a significant role in deter- 
mining galactic structure unless they either move galax- 
ies along the surface in {re,Ie,V) space or act in con- 
cert to preserve the manifold as the locus of equilibrium 
points of galactic structure. To reiterate, the important 
aspect of Figure 3 is not its linear nature, which is a re- 
sult of our assertion that a fundamental manifold exists 
for all galaxies, but rather the low scatter and lack of 
any systematic departures for specific classes of galaxies, 
both of which imply that the assumptions that we have 
made so far are appropriate to this level of precision. 
In other words, random or systematic variations of A 
and B across either the mass range of galaxies or galaxy 
types, variations from isothermality, and any other fac- 
tors that are not considered contribute scatter that is at 
most the observed scatter, which is comparable to that 
measured in TF or FP studies. We have achieved the 



first of the goals described in fJT] and now proceed to ex- 
amine whether our attribution of the departures from 
homology to Tg is correct. 

3.3. The Physical Validity of T{ 

The mathematical trick of placing all of the galaxy 
formation physics beyond the virial theorem into Tg po- 
tentially masks the importance of A and B. To check 
whether T{ truly reflects Tg, or whether it is in actu- 
ality a composite of various terms, we compare T{ to 
independent determinations of Tg. We do this for both 
normal ellipticals and dSphs. 

First, we compare to values of Te derived from a full 
Schwarschild analysis of the 2D line-of-si ght velocity dis- 
tribu t ions, T^'^'', of normal ellipticals (|Cappellari et al.l 
l2006l l2007t ). Because of the unknown constant in the 
definition of T{ relative to Tg, we have the freedom to 
normalize T{ to best match the ICappellari et al.l ()2006f ) 
data, which we do below (Figure |4]) . Figure |4] illustrates 
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Fig. 4. — A comparison of niass-to-light ratios derived from in- 
dependent means, X^^^ and our estimates of Te using Ti . For 
normal ellipticals (circles), we compare the mass- to-light ratio de- 
rived using the Schwarschild meth od of dynamical modeling for 
a set of nearby spheroidal galaxies l lCappellari et an i2006) to our 

estimates of Te using T{ . Open circles represent galaxies that 
[CapDcllari et al. (2006) note are fast rotators and filled circles rep- 
resent those that are not. The line is the 1:1 correspondence. The 
filled circles, which are the most appropriate comparison sample, 
show only 0.06 dex scatter (15% in T^). For the Milky Way dSph 
galaxies (triangles), we compare the mass-to -light ratio derived 
from fitting NFW profiles to kinematic data H Walker et al.ll2007l ) 

with our estimates of Te using Tl . The unknown constant relating 
Te to Te is set to ensure agreement in the mean values of X^^^ 
and Te and that value (—0.8) is then adopted for Eq. |H] 



the excellent correspondence between T{ and Tf'^^. The 
agreement is particularly good (0.06 dex rms, 15% in 
Te) for the galaxies that are most appropriate for our 
construction, namely those where the velocity disper- 
sion dominates over systemic rotation and anisotropy 
measures are small ( |/3, 7,(5| < 0T5 as measured by 
iCappellari et al.1 (120071 )1 . The scatter for those with large 
anisotropics is significantly greater (0.17 dex rms, 48% in 
Te), suggesting that a full knowledge of A and B would 
decrease the scatter in Figure [3] among those galaxies 
that do not fully satisfy the basic assumptions of our 
approach. 

Second, we examine whether this correspondence holds 
across the range of galaxy masses. In Figure H] we in- 
clude values of log Te for Galactic dSphs estimated us- 
ing NFW model fits to the extensive kinematic data of 
IWalker et all (|2007t l^. Unlike the Schwarschild analy- 
sis of normal ellipticals, which has the fr eedom to in- 
clude anisotropic velocity distributions, the I Walker et alj 
(J2007) analysis does not. The dSph data suggest a 
slightly different constant offset between logT{ and 
logTe (—0.78 rather than —0.82), but this change is 
modest and has a nearly undetectable effect on the scal- 
ing relation when ignored (see Figure [5|) . We adopt the 
average offset, —0.8, as the normalization constant and 
suggest an uncertainty in this number of ~ a few hun- 
dredths. Using this correspondence, we replace the con- 

** The calculated Te's use masses enclosed within r^ as calculated 
from the published fits, courtesy of Matthew Walker. 



Fig. 5. — The scaling relationship, Eq. \E\ using Xe iii place of 
Te. Data and symbols are as in Figure |4] The left panel includes 
all of the spheroids with T^^^, and the scatter is 0.09 dex about 
the 1:1 line with a mean offset of 0.005. In the right panel, we have 
removed the e llipticals with either significant rotation or anisotropy 
l|CaDpellari et al.. .2006 . .2007f ). and the scatter drops to 0.04 dex 
about the 1:1 liiijs with a ijiean offset of 0.004. 

stant m tq. [5jto obtain 

log re - log V^ + log /e + log Te + 0.8 = 0, (6) 

where one can either evaluate Te in some independent 
manner or express it in terms of V and le using the fit 
given in Table 1 and replacing logTg + 0.8 with logT{. 
The applicability of the same normalization constant 
for both normal ellipticals and dSphs supports the con- 
tention that structural variations, as reflected by changes 
in A and B, are modest over most of the mass scale 
covered in Figure [31 Using these independently derived 
measures of Te, we now return to Eq. [6l use the lit- 
erature values for Te rather than our fitting function, 
and plot the results in Figure O The difference between 
the left and right panels is the exclusion of the ellipti- 
cals that show evidence for rotation or anisotropic ve- 
locity dispersions. The data in both panels follow the 
1:1 correspondence well, although the correspondence in 
the right panel is striking. The scatter in that panel is 
0.04 dex, or less than 10% in the parameter values them- 
selves. We conclude that to a level of precision between 
10 and 25% (the scatter measured using these indepen- 
dently measured values of Te and the scatter measured 
using our fitting function, T{, respectively), Te encom- 
passes all of the additional physics necessary to proceed 
from the virial theorem to a description of galactic struc- 
ture. Thus, we achieve the second goal listed in fJTJ 

3.4. Evolving onto the Manifold 

For various reasons, the small scatter seen in Figures 
[3] and [5] is remarkable. Even if re and V are the same in 
two similar galaxies, one might expect variations in the 
stellar mass-to-light ratios, T*, of more than 50%, which 
would introduce scatter via variations in le- We suspect 
that at least part of the reason for the low observed scat- 
ter lies in the selection of galaxies in TF and FP stud- 
ies, which generally favor evolved, dynamically-relaxed 
galaxies, which are unlikely to exhibit dramatic varia- 
tions in T*. Therefore, we now return to the last of our 
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Fig. 6. — Results for the lGeha et al.l l|2006l ') sample of disk galax- 
ies with high gas mass fractions. Upper panel shows the galaxies 
(points) in comparison to the fitted relationship (line). Lower panel 
shows to where the galaxies might evolve as their gas is turned into 

stars. 



disk galaxy samples fGeha et al.ll200(t ). which contains 
disks with extremely high gas mass fractions, and might 
therefore be expected to harbor systems with dramati- 
cally different values of T* and Tg than those included 
in our analysis so far. 

When the majority of the baryons in a galaxy are in the 
gas rather than in stars, one might expect that our treat- 
ment as described above will fail because the connection 
between optical luminosity and mass via Te become s ten- 
uous. In fact, the galaxies in the lGeha et al.l ()2006D sam- 
ple do fall off the surface, as shown in the upper panel of 
Figure [6l However, a s demonstrated with regards to th e 
TuUy-Fisher relation (|McGaughll20"ol IGeha et al.l[200l . 
gas-rich and gas-poor galaxies have consistent scaling re- 
lationships if one considers total baryonic content instead 
of just that in the stellar component. Reviewing Eqs. [2] 
and [21 it is evident that the derivation of later equa- 
tions, such as Eq. [SI depended on a proxy for the en- 
closed mass. In light of previous studies, while optical 
luminosity may be an adequate proxy when most of the 
baryons are in stars, it is clearly inadequate when the 
majority of baryons are in the gas (or, more precisely, 
our fitting formula for Tg fails when different fractions of 
the baryons are luminous in otherwise similar galaxies). 
The reason optical luminosity works well for most disks 
is that the majority of their cold gas has been turned 
into stars (, Robcrts 1969; Bothun 1984; Kannappan 2004; 
[Read &:"^enthaml |2005[ ). Given that we do not have 
measured colors and gas masses for most of our sample, 
we cannot reformulate everything in terms of baryonic 
mass, but we can ask where these gas-rich galaxies would 
lie with respect to the 1:1 relationship in Figure [5| if they 
turned their gas into stars. Can they evolve onto the 
manifold? 

To complete this exercise, we make two questionable 
assumptions. First, we assume that r^. does not change 
during the conversion of gas to stars. Because the 
gaseous and stellar radial distributions are likely to be 
different, this assumption must fail at some level. Sec- 
ond, we assume that the stellar population formed from 



the gas eventually has the same stellar mass-to-light ra- 
tio (i.e., for stars only) as the late-type galaxies in our 
sample (T* = 0.97 in the J-band f or a typ ical Sbc galaxy 
with B~V ^ 0.57; iFukugita et a l. 1995; Bell fc de Joag 
120011 ). For these two assumptions, we then find the frac- 
tion of gas turned into stars that produces the best match 
to our scaling relation (Eq. [S]). We exclude one galaxy 
with a rotation velocity of 5 km sec~^, which appears to 
be unphysically small, and one galaxy with b/a > 0.9, 
for which it is difficult to deproject the velocity width. 

The result of this exercise, that these high gas mass 
fraction galaxies can evolve onto the relationship defined 
by the larger sample, is shown in the lower panel of Fig- 
ure[ni The resulting scatter is 0.16 dex, which is slightly 
lower than that found for the dSphs, but larger than for 
all of the galaxies combined. The difficulty with our sce- 
nario is that the agreement requires turning about 20% 
more gas than is available throughout the galaxy into 
stars. Perhaps this failure reflects the need for infalling 
gas, but it may also mean that as the galaxy evolves 
there is a corresponding change in Vg. For example, in a 
model where the gas is funneled efficiently into the center 
so that Te is only one-third of its current value, and 75% 
of the current gas is turned into stars, the scatter about 
the relationship is only 0.14 dex. We conclude that there 
are plausible scenarios in which these gas-rich galaxies 
lie on the observed relationship and that one might iden- 
tify galaxies that are still strongly evolving or those with 
a significant reservoir of cold gas as outliers from the 
manifold. Nevertheless, if these galaxies either convert 
their gas to stars or if we properly account for their en- 
tire baryonic content within r^, then we expect that they 
will satisfy the same scaling relation as all other galaxies. 

3.5. Revisiting our Simplifications 

Before we proceed to discuss further implications of 
these results, we step back to explore how a failure to 
satisfy the assumptions invoked in our simplifications of 
the virial theorem would manifest itself in our evaluation 
of Eq. [6] The potential "failures" fall into three classes. 

First, we might have introduced errors that are con- 
stant across the galaxy population. An example of such 
an error would be if we always underestimated the po- 
tential energy in our evaluation of the virial theorem by 
a fixed factor. Such an error would manifest itself as a 
zero point shift of the data relative to the expectation. 
Because we do not calculate the specific constant in Eq. 
[5] from any physical argument, this type of error will be 
transparently absorbed into the constant term when we 
determine it using independent measurements, as done 
in H3.3I For almost all of our discussion, this type of 
error is difhcult to detect but irrelevant. 

Second, we might have introduced errors that vary sys- 
tematically across the galaxy population. An example of 
such an error would be if we underestimated the poten- 
tial energy by a certain factor for low mass systems but 
overestimated it by a similar factor for high mass sys- 
tems. Such errors, to the degree that they correlate with 
at least one of the structural parameters, will lead to 
changes in the coefficients in Eq. [Sjor that describe T{ 
(Table [1]), but would not introduce scatter. This effect is 
analogous to introducing a "tilt" in FP analyses. Iden- 
tifying this type of error is critical if one aims to under- 
stand the specific nature of the fitted relationships such 
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as that describing T{ or to compare with simulations. 
We implicitly tested for such effects across galaxy types 
in ^3.21 and across galaxy mass in Wd.'d[ The coefficients 
derived from the virial theorem and a mass-to-light ratio 
that scales directly with independently-derived mass-to- 
light ratios successfully produce a tight scaling relation 
(Eq. [5] and Fig. [3|) . This result demonstrates that there 
is no effective "tilt" either with morphological type or 
across the full range of galaxy luminosities. 

Third, we might have introduced errors that are vari- 
able and not systematic across the galaxy population. 
An example of such an error would be if we have ignored 
a key determining factor of galactic structure, e.g. the 
number of nearby neighbors. In naive models, close pas- 
sages affect the luminosity of the system but do not af- 
fect the size or internal kinematics — leading to potential 
outHers in Eq. El Such effects, to the degree that they 
do not correlate with the remaining structural parame- 
ters, will introduce scatter at each point in the (rg, /g, V) 
space. The low scatter in both the entire sample (Fig. [3]) 
and the subsample with independently-measured mass- 
to-light ratios (Fig. [5]) demonstrates that any such errors 
introduce little noise. 

3.6. Connecting Tg to Physical Parameters 

What we have discussed so far are the end-products, 
or observables (rg, /g, V, Tg), rather than inputs, or 
true physical quantities, that determine the structure 
of galaxies. Two natural candidates for the parame- 
ters that drive galactic structure are the mass and an- 
gular momentum of a galaxy. For disks, analytic treat- 
ments of galaxy formation usi ng these two variables 
have been relatively s uccessful ([Fall fc EfstathioiilllQSOl : 
iDalcanton et al.iri997| ). Although these models require a 
few key assumptions that may not be entirely accurate, 
such as the conservation of angular momentum during 
collapse, their success, in combination with the results 
presented here, suggest that simple dynamical models 
may be able to reproduce the observable properties of 
galaxies. However, proceeding from such difficult to mea- 
sure quantities as mass and angular momentum to the 
observed structure of galaxies in a single step is likely to 
prove difficult. 

We focus instead on what we have learned so far from 
our analysis. Given that the virial theorem plus Tg are 
all one needs to generate a gross description of galactic 
structure fi j3.ll — 13. 3p . all of the interesting astrophysics 
of galaxy formation — at least as related to determining 
the current, gross structure of a galaxy — is encapsu- 
lated in Tg. What determines differences in Tg among 
galaxies? 

First, galaxies might convert a different fraction, ry, of 
their baryons to stars. Assuming a universal baryon-to- 
dark matter ratio for halos, this fraction can be mea- 
sured using the total mass-to-light ratio, T200, evalu- 
ated at i?200: the radius within which the mass den- 
sity is 200 times the critical density and the sys- 
tem is roughly virialized. Systematic variations in 77 
have already b een noted in studie s of halo occupa- 
tion statistics ()van den Bosch et al.l 120031: lY ang ct al.' 
2005f) ■ lensing mass meas urements ( Hoekstra e t al. 2005; 
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baum et al. 200 61), and direct bary on measure- 



ments (lLineLaL.20M; .GonzalezeLalJl2Q07l) . AU of these 
studies show that T200, and therefore r], depend on halo 
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Fig. 7. — Fraction of baryons that are converted to stars, rj, 
as a function of L/L* for spheroids (upper) and disks ( lower). 
Our d ata are the small dots. The values from Mandclbau m et al.l 
l|2006l) are shown as large filled circles with error bars and repre- 
sent averag e values for bins of L/L * . This plot includes the cluster 
spheroids of lZaritskv et al"] l|2006al ) as open triangles, so that we ex- 
tend the r ange of L/L* to the higher luminosities (see text) probed 
bv iMaiKlerbaum et al.l 1I2OOI) . 



mass. 

Second, the stars might be concentrated to varying de- 
grees relative to the dark matter due to differences in 
the assembly history. We quantify the stellar concentra- 
tion, ^, using ^ = i?20o/'"e. As the stars become more 
concentrated in the halo, ^ increases and Tg decreases. 

We now return to Eq. [2] with the aim of extracting 
from it expressions for the mass fraction of baryons that 
are converted to stars, 77, and the degree to which the 
stars are concentrated relative to the dark matter, ^. We 
rewrite Eq. [2] as appropriate at the virial radius, i?200i 



^200^200 



B200GM2 
R200 



(7) 



where M200 = 37i'i?20oP20o, P200 = 200pcrit, and R200 = 
^rg, where Pcrit is the universal critical mass density at 
the present epoch. To make further progress, we set 
V — V200. This is patently incorrect both because the 
dark matter potential itself is u nlikely to be isothermal 
out to i?200 (jNavarro et al.l[l997b and any central concen- 
tration of the baryons will affect V. However, any con- 
stant fractional differences — for example, V = I.2V200 
for all galaxies — will be absorbed later into our nor- 
malization. What do concern us, but are ignored here, 
are differences in this veloc ity scaling that depe nd on 
the properties of the galaxy (jCourteau et al.|[2007[ ) . This 
problem might be correctable in an iterative scheme (i.e., 
assume a non-isothermal potential, estimate ^, evaluate 
the difference between V and V200 , correct V200 , and iter- 
ate until convergenc e) or in a more sophist icated model 
of galaxy formation (jSomerville et al.l[2007| ). but both of 
these remedies require some model assumptions and lead 
us from the analytic descriptions that we aim to explore. 
Work is needed to determine the magnitude of the error 
introduced by our simple treatment. 
Continuing, we define rj through the equation 
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L = 



tiIbM, 



200 



(8) 



where T^,, is the mass-to-hght ratio of the stellar popu- 
lation, L is the total luminosity, and Jb is the baryon 
mass fraction. Algebra enables us to derive equations for 
^ and J]: 
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(10) 



We express the combination of these two quantities as 
log C + log 77 = log— hlog— j^, (11) 



^3 



where all constants, as well as the structural factors A200 
and -62001 in Eqs. [5] and (TU] are contained in K and 
J . ^200 E^nd i3200 are not necessarily constant, although 
their analogs at r^ are well-behaved f ij3.2p . We assume 
that j42oo and i?200 are similarly well-behaved at -R200- 
Rewriting, we get 

log^-|-log?7 = — 2l0gF-f logTe +log/e + log T, + COIlst. 

(12) 
We know from t ^3.2l that we can replace the three leading 
terms on the right hand side with — log Te to within a 
constant, so 

log ^ -I- log 1] = — log Te -I- log T* -|- const. (13) 

Because Te — f{V,Ie), as defined in Table 1, log?7 + log^ 
is also a function V and le- 

To provide a more direct example of the possi- 
ble use of these equatio ns, we use the results from 
iMandelbaum et aTl ()2006f ) to evaluate the constant terms 
(including the ass umed constant terrns Aan n and i320o) 
in Eqs. [H]and[Tni IMandelbaum et al.l ()2006l ) use results 
from weak lensing to evaluate the fraction of baryons 
that are turned into stars as a function of galaxy lumi- 
nosity. They provide empirical values for 77 as a function 
of galaxy luminosity and morphology. We use their re- 
sults for f] for L* spirals, and assume T, = 1.7 in the 
/-band and a unive rsal baryon mass fraction of 0.175 
(jSpergel et al.l [20071 ) . to set the values of the leading 
coefficients in our Eqs. [9] and [TOl ^ = lAVr~^ and 
77 = 1.9 X 10-5(L/L*)T*V^-3. The adopted value of T» 
for an L* spir al (assuming that it is an Sab type with 
B - V= 0.8; (Fukugita et al."1995')) is calculated using 
Table 1 of^BcU & de Jong (2001). We then plot our cal- 
culated values of -q for all of our galaxies (Figure [7|) , as- 
suming that T« is 1.7 and 2.5 for disks and spheroids, 
respectively. For the purpose of this comparison, we aug- 
ment our galaxy sample with measurements of cluster 
spheroids (CSph), the brightest cluster galaxy plus the 
intracluster stars of groups and clusters (jZaritskv et al.l 
[2006a) . Our previous work shows that these systems lie 
on the FM and, in this context, they allow us to extend 
the range of L/L* over w hich we can compare to the 
IMandelbaum et aD (|2006( l results. We decided against 
including the CSphs throughout our current study be- 
cause our focus is on galaxies, but CSphs do indeed fall 
on the 1:1 line in Figure [3l 



Our values for r/ are in agreement over the 
range of luminosit i es an d galaxy types presented by 
IMandelbaum et al] (|2006D . except possibly for the high- 
est luminosity spheroids: brightest cluster galaxies and 
CSphs. Our data suggest a qualitativ ely similar drop in 77 
to that found bv IMandelbaum et al.l (J2006) . although at 
somewhat larger values of L/L* . There are many tech- 
nical reasons (such as the use of different bandpasses, 
subtleties in the definition of total magnitudes, and com- 
plications introduced by intracluster light for these most 
massive spheroids) that preclude any conclusion about 
whether there is a true discrepancy. In general, we 
agree quite well both qualitatively and quantitatively 
with their results. This agreement, in turn, suggests that 
yl2oo and i?200 do not vary strongly as a function of mass 
or morphological type. 

We proceed now to calculate 77 and ^ as a function of 
V for all of our galaxies and show the results in Figure 
m Here, to within the fiaws in our simple derivation and 
the heterogeneous sample, is a full description of how 
baryons turn into stars and distribute themselves in all 
galaxies ranging from dSphs to BCGs. There are several 
striking results. First, the relatively large values of Te 
for dwarf spheroidals are primarily driven by rj rather 
than by ^, which is surprisingly constant across the full 
range of V . Second, at a given V , ellipticals and spirals 
have similar values of 7;. Therefore, the difference in their 
values of Te are due primarily to differences in ^. Third, 
there is a steady decline in 77 for logF > 1.5 {V > 32 km 
s~^) even for the spheroids, among which large variations 
in T, are less likely. Systems with comparable V to 
the Milky Way, 149 < 1^ < 163 km s~^ or alternatively 
210 < Vc < 230 km s~\ have 77 = 0.14 ± 0.04, with the 
spheroids being on average 5 times more concentrated 
than the disks. All of these results await resolution of 
two key open questions in the evaluation of ^ and 77: 1) 
how to treat the difference between V and V200 and 2) 
whether the structural terms are as well-behaved at R200 
as they are at r^. Nevertheless, this analysis illustrates 
how we might construct a bridge between the empirical 
relations based on observables and more theoretical ones 
based on fundamental physical parameters. Thus, we 
achieve the third goal listed in SJTJ 

4. SUMMARY 

We have shown that all classes of galaxies, ranging in 
mass from dwarf spheroidals to brightest cluster galax- 
ies, and in type from spheroids to disks, fall on a two 
dimensional surface in the observable space {ve, Ie,V), 
where V^ = ^v^ + a^ , over three orders of magnitude in 
Te with < 25% scatter. The scatter about that surface is 
comparable to that observed in Fundamental Plane and 
TuUy-Fisher studies in which the range of galaxy types 
and luminosities is much more limited. The TF and FP 
relationships are subsets of the manifold presented here. 
This finding alone demonstrates that the structure of all 
galaxies can be described with a highly limited set of pa- 
rameters. The observational ones, (re,/e,^), may not 
be optimal, even though they do a remarkably good job. 
The small scatter about the mean relation implies that 
a host of potential physical phenomena such as environ- 
mental effects, star formation history, nuclear activity, 
accretion history, and feedback are either (1) relatively 
unimportant in determining the structure of galaxies, (2) 
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Fig. 8. — Mass fraction of baryons that are converted to stars, 77, (top) and stellar concentration, 5, (bottom) as a a function of V for 
our entire sample (spheroids represented with filled circles, disks with open circles). Various results are in evidence, including the dramatic 
decrease in rj for the dSphs, the similarity in ri among spirals and spheroids in the regime where they overlap in V , the systemic decline in 
r) with increasing V for logV > 1.5, and the generally greater § for spheroids relative to disks in that same velocity range. 

move galaxies along this well-defined relationship, or (3) 
balance each other so as to define the mean relation as 
the locus of galactic structure equilibria. 

We developed a simple analytic treatment in which 
we asserted the existence of a fundamental manifold of 
galaxies. By requiring the simple virial theorem deriva- 
tion to result in a two dimensional manifold in observed 
space, we specify the behavior of the mass-to-light ratio 
within Te, Tg. We then tested this assertion by compar- 
ing our inferred values of Te to those derived indepen- 
dently from much more sophisticated modeling for both 
normal ellipticals and dSphs. The agreement is quanti- 
tatively excellent, with less than 15% scatter in mass-to- 
light ratios for those galaxies that satisfy our dynami- 
cal criteria. This result demonstrates that the principal 
additional ingredient necessary in proceeding from the 
virial theorem to a description of galactic structure is 
knowledge of the mass-to- light ratio within re. Addi- 
tional factors, which could have been important, such as 



internal kinematic anisotropy or variations in the radial 
profile of the gravitational potential from one galaxy to 
another, must play a role at less than the 25% level. The 
observed manifold is described by 

log re - log y2 + log l^ + log Te + 0.8 = 0, (14) 

where we also provide a fitting function for logTe in 
terms of V and le- The equations presented here are 
numerically appropriate for re, V ^ /e, and T in units of 
kpc, km s~^, Lo/pc^^, and solar units, respectively, and 
based primarily on /-band observations. 

We then discuss what the inferred behavior of Tg may 
mean for the physical characteristics of the galaxies. In 
particular, we speculate that the two principal determi- 
nants of Te arc the mass fraction of baryons that are 
turned into stars, rj, and the degree to which the stars 
are spatially concentrated relative to the dark matter, 
^ = i?20o/^e- We derive equations for both quantities in 
terms of unknown structural parameters and the observ- 
ables. We relate the two quantities using the expression 
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that we de rived for Te- Finally, we u se independent mea- 
sures of rj (jMandelbaum et al]l2006( l to solve for the un- 
known structural terms for one set of galaxies and then 
compare the behavior of rj across other luminosities and 
galaxy types as determined both from our analysis and 
that independent weak lensing study. This comparison 
leads to simple expressions for r/ and ^, 



7? = 1.9x 10 



_5(i/i*)T, 



y3 



and 






(15) 
(16) 



We are then able to extend the measurements of ry and 
^ to the full range of galaxies. As rough guides, we find 
that, for most galaxies, 0.04 < 77 < 0.6 and 10 < ^ < 200, 
although these can be evaluated on a galaxy-by-galaxy 
basis. Systems with comparable V to the Milky Way, 
149 < T^ < 163 km s'^ or alternatively 210 < Vc < 230 



km 



have 77 = 0.14 ± 0.04, with the spheroids being 



on average 5 times more concentrated than the disks. 
Overall, we reach a set of general conclusions. First, the 
relatively large values of Te for dwarf sphcroidals are pri- 
marily driven by 77 rather than by ^, which is surprisingly 
constant across the full range of V. Second, at a given 
V, ellipticals and spirals have similar values of 77 (< 10% 
difference for spheroids and disks with 149 < V^ < 163 
km s~^). Therefore, the difference in their values of Te is 
due primarily to differences in ^. Third, there is a steady 
decline in 7; for logF > 1.5 (1^ > 32 km s^^). 

The data used here fall short of the ideal sample 
from which to properly derive the quantitative values 
that mathematically describe the manifold, primarily 
due to the heterogeneous nature of the amalgamated 
sample. Nevertheless, the sample does demonstrate that 
the range of galaxy structure is dominated by only two 
parameters. The lower scatt er obtained eith e r for a sin- 
gle sample (0.06 dex for the ISpringob et ahl (2007) disk 
sample) or for i ndependently derived Te 's (0.04 dex when 
using both the ICappellari et all ()2006f ) data for ellipti- 
cals and the I Walker et al.l (|2007| ) data for dSphs) suggest 
that a homogeneous sample might show that the myriad 
of possible influences on galactic structure (environment, 
accretion history, AGN activity, star formation history, 
and others) contribute at most a ~ 10% scatter to the 
scaling relationship presented in Eq. [Ml 



The existence of a highly constrained surface on which 
galaxies lie does not eliminate the need for additional 
physics. In particular, as we have hinted throughout, 
many physical effects might move galaxies along the sur- 
face or perhaps counter-balance to move galaxies back 
to the equilibrium surface described by the manifold. As 
such, future galaxy models may be more constrained by 
the distribution of galaxies on the surface rather than 
perpendicular to it. Our heterogeneous sample is ill- 
suited to say much about the distribution of sources on 
the surface. Much work still remains. 

We close by returning to the analogy of stellar struc- 
ture. It is evident that we are still far from a physical 
theory of galactic structure, but that we have progressed 
in several key aspects. First, we have now demonstrated 
that the entire family of galaxies can be described by 
sets of two parameters (e.g., V and le or 77 and ^). This 
finding motivates the search for relatively simple expres- 
sions of galactic structure that are connected to a small 
set of physical parameters, such as mass and angular mo- 
mentum. Second, we have identified the principal char- 
acteristic that remains to be explained, namely Te. The 
virial theorem plus an understanding of Te are all that 
are necessary to predict the size, internal kinematics, or 
luminosity of a galaxy, when given the other two. This in 
turn places the focus on understanding what determines 
the fraction of baryons that are turned into stars and how 
those stars are packed within the dark halo. If those 
quantities can then be connected to more fundamental 
parameters, such as mass and angular momentum, then 
one could proceed from the physical parameters directly 
to the observables. At that point, we will have indeed 
produced equations of galactic structure. 
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